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Abstract. We analyze a 137 ks exposure X-ray spectrum of the Seyfert 1 galaxy NGC 5548 obtained with the XMM-Newton 
Reflection Grating Spectrometer. Due to the long exposure time, the spectrum is of higher statistical quality than the previous 
observations of this AGN. Therefore, we detect for the first time in NGC 5548 inner-shell transitions from O III to O VI ions, and 
the Unresolved Transition Array of M-shell iron. The warm absorber found from this X-ray observation spans three orders of 
magnitude in ionization parameter. We detect O III, which is as lowly ionized as the warm absorber detected in the UV band, to 
Fe XXIV. For O VI the column density determined from our X-ray data is an order of magnitude larger than the column density 
measured in previous UV observations. We conclude that there is substantially more low ionized material than previously 
deduced from UV observations. However, only a few percent of the warm absorber detected in the X-rays is lowly ionized. 
A 99.9% significant increase in the derived absorbing column density with higher ionization states is observed. The outflow 
velocities determined from the X-ray absorption lines are consistent with those deduced from the UV lines, evidence, together 
with the detection of O VI, that the X-ray and UV warm absorber are different manifestations of the same phenomenon. From 
a simple mass conservation argument, we indicate that our data set is consistent with an outflow with small opening angle 
formed due to instabilities in the accretion disk. Possible due to uncertainties in the radiative transport mechanism, an apparent 
deviant iron to oxygen abundance is detected. No strong relativistically broadened emission lines of O VIII, N VII and C VI were 
detected. 
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1. Introduction 

The immediate environment of active galactic nuclei (AGN) 
is poorly understood. The black hole is fed from an accre- 
tion disk, which becomes visible through high energy radiation. 
Possibly powered by radiation pressure, an outflowing wind is 
formed from the accretion disk. In some AGN, the wind is vis- 
ible through its continuum and line absorption of the radiation 
from the nucleus. These absorption lines were already detected 
in the UV part of the spectrum (Mathur, Elvis & Wilkes 1995 1, 
and have now been confirmed in the X-ray spectra (Kaastra et 
al. l2000l i. 

NGC 5548, a Seyfert 1 galaxy, has been well studied in the 
X-rays due to its relative X-ray brightness, and the fact that it 
is a nearby active galactic nucleus with a redshift of z=0.01676 



(Crenshaw & Rraemer 1999). This Seyfert galaxy was stud- 
ied in detail with ROSAT (Done et al. U995> . They detected an 
O VII and O VIII K-shell absorption edge, and concluded that 
there was a warm absorber present. Therefore, NGC 5548 is an 
ideal Seyfert 1 galaxy to study with the high resolution spec- 
troscopy instruments on-board of XMM-Newton and Chandra. 
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Tabled list the observations made with these instruments. In 
this paper we discuss a 137 ks exposure time Reflection Grating 
Spectrometer (RGS) spectrum from the guest observation pro- 
gram of XMM-Newton. 

2. Observations and data reduction 

Here we present data taken with the RGS instrument (den 
Herder et al. 2001 ). The EPIC data (the continuum spectrum, 
and the Fe Ka line) are presented separately (Pounds et al. 
2002 ). The time variability of the combined XMM-Newton and 
RXTE data sets will be presented by Markowitz et al. (2002). 

Shortly after launch two of the eighteen CCD chains failed. 
For RGS 1 this causes a data gap in the 10 -14 A range (where 
most Fe-L and Ne absorption lines are), for RGS 2 the gap is 
between 20 - 24 A (where the O VII triplet is). 

The data were reduced using the developer's version of the 
XMM-Newton data analysis package SAS version 5.3, which 
comprises a correct calibration of the effective area and the in- 
strumental O-edge. The latter is important as the ls-2p absorp- 
tion lines of O II and O III are cosmologically redshifted toward 
the instrumental O I edge. 

In this paper we analyze the first and second order spectra 
obtained from both RGS 1 and 2. Although the second order 
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Table 1. High resolution X-ray observations of NGC 5548. 
LETGS: Low Energy Transmission Grating in combination 
with the High Resolution Camera; HETGS: High Energy 
Transmission Grating in combination with the ACIS camera 
both onboard Chandra; RGS: Reflection Grating Spectrometer 
onboard XMM-Newton. The exposure time is given in in ks. 



Date 


Instr. 


Expos. 


References-comments 


1999 Dec. 11 


LETGS 


86.4 


Kaastra et al. 2000, 2002 


2000 Feb. 5 


HETGS 


83 


Yaqoobetal. 2001, 








Kaastra et al. 2002 


2000 Dec. 24 


RGS 


28 


noisy spectrum 








low flux state 


2001 Jul. 09-12 


RGS 


137 


present work 




10 15 20 25 30 35 

Wavelength (A) 



Fig. 1. The fluxed spectrum of NGC 5548, RGS 1 and 2 are 
added together and no error bars are plotted for clarity. The 
typical 1 g error on the data is 1 photon m~ 2 s _1 A -1 . Some 
of the strongest features are identified. 



spectrum has a lower count rate and only a wavelength cover- 
age from 5 to 19 A, it does have a higher wavelength resolution 
than the first order. Throughout the paper we have fitted RGS 1 
and 2 simultaneously, but for clarity of representation we have 
added the RGS 1 and RGS 2 spectra in the figures. 

All the spectral analysis was done with the SPEX pack- 
age (Kaastra et al. 2002bi. Fig.^shows the fluxed spectrum of 
NGC 5548, indicating the Unresolved Transition Array (UTA) 
of inner-shell iron absorption lines and some other prominent 
features. 

2. 1. Wavelength scale accuracy 

In general the RGS wavelength scale can be reconstructed to 
an accuracy of 8 mA (den Herder et al. 2001 1. However, in 
our case larger uncertainties exist, due to errors in the recorded 
satellite attitude and right ascension in the data files available. 
From the analysis of the first order data we find that the O VII 
forbidden line is blueshifted by 285 km s" 1 (-0.021 A). The 
O VII forbidden line is not blueshifted in the earlier LETGS 
data of NGC 5548 (Kaastra et al. I2002a> . No blueshift for the 



forbidden line can be determined in second order. However, 
comparing first and second order results for the warm absorber, 
we find a smaller outflow velocity in the second order data set. 
The difference in wavelength between first and second order, 
for the strongest line in the second order, O VIII Ly/3 (both or- 
ders were fitted with the xabs model, see Sect. 3) is —0.019 
± 0.005 A. Finally, in the first order data set there is a weak 
Galactic O VII recombination line, which has a rest wavelength 
of 21.6 A. This absorption line also shows an apparent veloc- 
ity shift of -230 km s" 1 (-0.017 A), although no blueshift 
should be expected. To correct for this wavelength calibration 
problem we applied a constant wavelength shift of —0.017 A to 
the model for the first order spectra. If not corrected, this wave- 
length shift would result in a velocity error between 130 km 
s _1 for the longest wavelengths to 620 km s _1 for the shortest 
wavelengths. The residual error in the wavelength scale should 
be within the 8 mA uncertainty of the calibration. All velocities 
quoted in this paper are corrected for this wavelength shift. 

3. Modeling the warm absorber 

In NGC 5548 there is an intricate pattern of densely spaced 
absorption lines as well as the possible presence of weak rela- 
tivistic emission lines (Kaastra et al. 2002a I. Blending, mostly 
of iron lines, is important in many wavelength ranges, making 
equivalent width (EW) estimations rather complicated. On the 
other hand, in the X-ray band many of the observed ions have 
more than one detected absorption line, significantly improving 
the column density and velocity outflow measurements. This is 
also important in order to disentangle possible partial cover- 
ing factors and column densities of the ions. The possibility of 
relativistic emission lines complicates the analysis of the con- 
tinuum spectral shape, which is already complicated to derive 
due to the presence of many edges and absorption lines. We 
consider here two approaches for fitting the warm absorber. In 
the first approach all absorption lines and edges of each indi- 
vidual ion are calculated, using a separate value for the column 
density for each ion. The transmission of all ions is then com- 
bined to calculate the total transmission of the warm absorber. 
Outflow velocity and velocity broadening are coupled for the 
different ions. This is the slab model in SPEX (Kaastra et al. 
I2002b> . 

In the second approach the column densities of the ions 
are coupled through a photoionization model. In this case only 
the total column density, ionization parameter, elemental abun- 
dances and the outflow and broadening are free parameters. 
This is the xabs model. For a full description of both methods 
see Kaastra et al. ( 2002b I. 

The results obtained with the xabs model are more depen- 
dent on theory, because the ionization parameter determines 
the column density ratios of the ions. An advantage of the 
xabs model is that different ionization components, each with 
their own velocity structure can be combined. This is harder 
to accomplish with slab components, in particular if the veloc- 
ity components are barely separable with the current spectral 
resolution. This is the case in NGC 5548, and would lead to 
strongly correlated parameters in the spectral fitting procedure. 
One further difference between the models is that ions which 
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Table 2. The parameter values quoted are those obtained in the 
xabs model fitting and differ from the best fit continuum with- 
out absorption. The continuum model was used as a basis for 
adding the warm absorber, either with a slab component (see 
Table [3} or with xabs components (see Table [4). The quoted 
Xv 2 is for the best fit continuum parameters without absorp- 
tion, and thus not the parameters quoted in the table, which we 
derived including absorption. 



X 


j 2 equals 


1.96 for 2441 degrees of freedom 


pi: 


norm" 


(1.57 ± 0.03) x 10 52 ph s" 1 keV" 1 




r 


1.77 ±0.02 


mbb: 


norm 6 


(6.5 ± 0.7) x 10 32 m 1/2 




T 


(97 ± 6) eV 



a at 1 keV. 



6 norm=emitting area times square root of electron density. 



have a rather small cross-section are automatically included in 
xabs, but since these are more difficult to discern individually, 
they are not usually included in the slab modeling. 

4. Data analysis 

4.1. The continuum 

The continuum model consists of a power-law spectrum with 
Galactic absorption, and a modified black body (see Kaastra 
& Barr 1989 for more details). A second modified black body 
with a high temperature is required to describe the EPIC data, 
but is not detected in the limited RGS band (Pounds et al. 
2002 1. This indicates that RGS data alone cannot constrain the 
continuum accurately. Table[2]gives the continuum parameters 
used for the analysis of the RGS data. 

The temperature (T) of the modified black body was al- 
lowed to vary as well as the normalizations of the power-law 
and the modified black body. However, the photon index (r), 
which is very sensitive to small changes in the model, as well 
as the Galactic H I column density (1.65 x 10 24 m~ 2 , Nandra 
et al. 119931 1 were kept fixed throughout the fitting process. 
The obtained parameters are within calibration uncertainties 
consistent with those obtained from simultaneous EPIC and 
BeppoSAX observations. 

4.2. Spectral fit using the slab model 

As a first approach to the analysis of the warm absorber, we 
modeled it by adding a slab component to the continuum 
model. We refitted the data, leaving the same parameters free 
as above for the continuum, as well as the overall outflow ve- 
locity, and all the relevant ion column densities. Finally, after 
obtaining this overall fit we also let the photon index vary, and 
now obtain a reasonable result. The results of this fit were then 
used to fine-tune the column density and outflow velocity of 
each ion which has a detectable absorption line or edge in the 
RGS wavelength range. To do so we decoupled subsequently 
for each ion its outflow velocity from the outflow velocity of all 
the other ions in the model. This method was followed to mini- 
mize the effects from blending. The significance of the different 



ionic columns can be derived from the best fit column densities 
and their error bars. The result of this fitting is summarized in 
Table |3] This Table also lists the ionization state for the maxi- 
mum column density of the particular ion. Note however, that 
each ion exists for a range of ionization states, and thus the 
quoted ionization is only an indication. Adding the slab com- 
ponent significantly improves the Xv 2 from 1.96 to 1.26 (with 
2404 degrees of freedom). 

For the ions with only an upper limit for their column den- 
sities, absorption is only detected through continuum absorp- 
tion (Ar XIV, Fe V and Fe VII) or a single/few weak lines (O IV, 
Si X, Si XI, Fe I, Fe XII and Fe XIV). As a result the column den- 
sity as well as the velocity of these ions are poorly determined. 
The other ions with poorly determined velocity shifts, but with 
solid column densities, are those for which the column density 
is determined from an edge (C IV, N V and Si IX), for which the 
absorption lines are strongly blended (Fe XII), or because only 
a few weak absorption lines are detected. 

We find a very well constrained column density for C IV 
and N V, although both these ions have in the RGS band only 
an absorption edge, at 36 A and 25 A (rest wavelength) respec- 
tively. As the depth of the edge correlates with the assumed 
underlying continuum, these column densities are rather more 
uncertain than the fit results indicate. To ascertain how these 
column densities change the overall continuum shape we refit- 
ted the data, but now forcing a column density for C IV and N V 
of zero. For the new fit, which gives a similarly good x^ 2 , the 
normalization of the modified black body increased by an order 
of magnitude, while the temperature of the same component 
decreased by an order of magnitude. As a result of this change 
in the continuum model, the column densities found for O III 
and Ne IX also changed significantly. In this respect the xabs 
modeling is more reliable, as the column densities determined 
for ions that are only detected through an edge, are based upon 
the determined ionization state. 

The same method for modeling the warm absorber was 
used for fitting the second order data. The second order data for 
RGS 1 and 2 were not fitted together with the first order data, 
as the systematic shift in the first order spectra is corrected for 
in the applied model and not in the data itself. However, both 
orders are in excellent agreement, see Fig. [2] 

The average outflow velocity found, including only those 
ions with well determined velocities in both orders, is —350 
± 250 km s _1 for the first order data set, in comparison with 
—420 ± 340 km s _1 for the second order data set. For fur- 
ther analysis we focus on the first order data, because the better 
statistics outweigh the better spectral resolution of the second 
order data. 

4.3. Spectral fit using the xabs model 

Using the same continuum model as in the previous section, 
we replaced the slab model by the xabs model for the warm 
absorber. First we fitted only the ionization parameter (£), the 
hydrogen column density {Nn) and the outflow velocity (v), 
in addition to the continuum parameters that we leave free 
throughout the fitting procedure. The abundances were kept to 
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Table 3. Best fit values for the logarithms of the column densities in mT 2 and the outflow velocity, using the slab component. 
The logarithm of the ionization parameter for which the ion has the largest column density as determined from XSTAR is given 
in column 4. 





log N H 


V 


log? 




log N H 


V 


log? 


ion 


(rrr 2 ) 


(kms^ 1 ) 


(10" 9 Wm) 


ion 


Cm" 2 ) 


(km s _1 ) 


(10 -9 W m) 


ClV a 


21.4 ±0.1 





-0.6 


Sxiv 


20.1 ± 0.3 


+90 ± 280 


2.4 


Cv 


20.4 ± 0.2 


-490 ± 200 


0.2 


ArXIII 


19.4 ± 0.4 


-630 ± 210 


1.8 


Cvi 


21.1 ± 0.1 


-510 ±70 


1.2 


ArXIV 


< 19.4 


_ 


2.1 


Nv° 


20.6 ± 0.3 


_ 


0.0 


Arxv 


20.1 ±0.3 


-470 ± 400 


2.3 


Nvi 


20.1 ± 0.1 


-470 ± 130 


0.7 


Fel 


< 19.8 


_ 


< -4.0 


N VII 


20.5 ±0.1 


-390 ± 150 


1.5 


Fe II 


19.9 ±0.3 





- 3.0 


OlII 


20.0 ± 0.6 





- 1.8 


Fev 


< 19.7 


_ 


- 1.3 


OlV 


< 20.2 


_ 


-0.7 


Fe VI 


19.8 ±0.4 





-0.9 


Ov 6 


20.7 ± 0.3 


+430 


0.0 


FeVII 


< 19.7 




-0.3 


Ovi 


20.5 ±0.1 


-250 ± 160 


0.5 


Fe VIII 


20.1 ±0.1 


-340 ± 450 


0.1 


Ovn 


21.57 ± 0.05 


-450 ± 90 


1.2 


FelX 


20.1 ±0.1 


+ 130 ± 220 


0.4 


O via 


22.24 ± 0.03 


-310 ±70 


1.8 


Fex 


19.9 ±0.1 


-670 ± 430 


0.6 


NelX 


21.82 ± 0.06 


-30 ± 200 


1.7 


FeXI 


20.1 ±0.1 


-430 ± 400 


0.8 


Nex 


21.5 ± 0.2 


+10 ± 230 


2.3 


Fe XII 


< 19.7 




1.0 


Nax 


20.5 ± 0.4 




1.9 


FeXIII 


19.4 ± 0.6 


_ 


1.2 


Naxi 


20.6 ± 0.6 




2.4 


FeXIV 


< 19.5 


_ 


1.4 


Mgx 


21.4 ±0.1 


-560 ± 250 


1.5 


FeXV 


< 19.7 


— 


1.6 


MgXI 


20.5 ± 0.5 




2.1 


FeXVI 


19.7 ±0.3 


-500 ± 450 


1.6 


MgXII 


20.8 ± 0.5 


-970 ± 990 


2.6 


FeXVII 


20.1 ±0.1 


-480 ± 170 


2.1 


Si IX 


21.3 ± 0.1 




0.9 


Fe XVIII 


20.1 ±0.1 


-420 ± 240 


2.3 


Six 


< 20.4 




1.3 


FeXIX 


20.49 ± 0.09 


-350 ± 230 


2.5 


Si XI 


<20.1 




1.7 


FeXX 


20.4 ± 0.2 


-960 ± 710 


2.8 


Si XII 


<20.5 


-400 ± 380 


2.1 


Fe XXI 


20.5 ± 0.2 


-65 ± 550 


3.0 


Sx 


20.3 ± 0.2 


-820 ± 420 


1.1 


Fe XXII 


20.6 ± 0.2 




3.2 


SXII 


20.3 ± 0.1 


-410 ± 80 


1.7 


Fe xxm 


20.4 ± 0.4 




3.3 


Sxin 


20.1 ± 0.2 


-640 ± 200 


2.1 


Fe XXIV 


20.7 ± 0.5 




3.5 



a See text for details. 

6 Possible error in the wavelength. 



produces no X-ray lines, we normalized the abundance rela- 
tive to oxygen, which was frozen to solar abundance. Finally 
we added two more xabs components to the model, both with 
lower ionization parameters, significantly improving the spec- 
tral fit. 

Adding the first xabs component improved \ v 2 from 1.97 
to 1 .66. The addition of the second xabs component, which was 
the one with the lowest ionization parameter, further improved 
Xv 2 to a value of 1.46. However, for this xabs component we 
found two abundances which have unphysical values. Namely, 
Ne and S are overabundant by a factor of 43 and 106, respec- 
tively. From detailed fitting, it is clear that both these abun- 
dances are mainly determined from the absorption edges in 
the spectrum, and are thus rather uncertain. Fitting the model 
again, but fixing the abundance of Ne and S in the lowest ion- 
ization component to solar, we find a similar good fit, namely 
with Xv 2 of 1 .48. The addition of the last xabs component gives 
a final Xv 2 = 1-22. Table|4]lists the best fit parameters, with the 
xabs components labeled to be consistent with the components 
found by Kaastra et al. ( 2002a I. The significance, according to 
F-tests, is at least 100, 99.8 and 99.9 percent for components A, 
B and C, respectively. Only in the last step we left the photon 



Table 4. The best fit results for the first order data xabs fit. 
Further details are given in the text. 



Component 


A 


B 


C 


logC(10"' J Wm) 


2.69±0.04 


1.98±0.06 


0.40±0.03 


log Nh (m -2 ) 


25.68±0.05 


25.52±0.06 


24.15±0.01 


outflow v (km s _1 ) 


-311±60 


_440±100 


-290±70 


abundances" : 








C 


<4.9 


0.6±0.2 


1.8±0.6 


N 


2.2±1.0 


<0.3 


1.7±0.5 


Ne 


2.5±1.4 


0.1 ±0.7 


b 


Na 


<8.3 


3.9±3.9 


<350 


Mg 


1.0±0.8 


<0.6 


<5.8 


Si 


2.2±2.2 


<0.2 


<3.9 


S 


0.6±0.7 


0.5±0.2 


b 


Ar 


2.2±2.2 


0.4±0.5 


<5.4 


Fe 


0.70±0.06 


0.11 ±0.03 


7.4±0.6 



a The abundances are relative to O, which is kept at solar value. 
b Frozen to solar value, see text for explanation. 



the solar values (Anders & Grevesse 1989). Later, the abun- 
dances were left as free parameters. However, as hydrogen 
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13 14 15 16 

Wavelength (A) 



Fig. 2. Best xabs fit to the first (top) and second order (bottom) 
data. First order data are rebinned by a factor of 2, second order 
by a factor of 3. For the first order data a CCD gap related 
feature can be seen at 13 A and a dead pixel is seen short- ward 
of 16.5 A. Between 12.5 and 14 A only RGS 2 gives results, 
hence the larger error bars. For second order data a dead pixel 
is seen at 14.05 A and a CCD gap related feature can be seen 
at 14.45 A. 

Table 5. Values determined for the forbidden line of O VII. 



Vllf: 


EW 


(-170±22)mA 




flux 


(0.9±0.1)phm -2 s _1 




FWHM 


< 36 mA 






<490kms~ 1 




A a 


(22.096 ± 0.005) A 




v b 


(-110 ±80) km s -1 



a Wavelengths are given in the rest frame of NGC 5548. 

6 Velocity shift from comparing rest and observed wavelengths. 



index free, but this did not improve the \v 2 ■ The continuum 
parameters for this model correspond within 1 a with the con- 
tinuum parameters as found from the slab fit. 

4.4. Emission and absorption by oxygen 

Fig EJshows the spectrum around the O VII triplet. In Table[5] 
we list the parameters of the O VII forbidden line in emission. 
For the Ovil forbidden line Table [6] gives the fluxes, outflow 
velocity as well as an upper limit to the width of the line, as 
measured from the HETGS and LETGS spectra (data taken 
from Kaastra et al. 2002a I and our data set. Note that this line 
has the same measured flux over a period of 1 .5 years, despite 
large variations in the continuum. In all three observations the 
line is unresolved. The above indicate that the O VII forbidden 
line is formed at least 0.46 pc from the ionizing source. 

Interestingly, the inter-combination line is not present in 
the data. An excess is observed at 21.74 A rest wavelength, 
shifted by 0.06 A from the rest wavelength of the inter- 
combination line. The relative wavelength accuracy of RGS is 
2 mA. Assuming that the inter-combination and forbidden lines 
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Table 6. Comparison of the O VII forbidden line from three 
different observations. 



Instrument 


V 




flux 




(kmr 1 ) 


(km s" 1 ) 


(ph m~ 2 s" 1 ) 


HETGS 


-70 ± 100 


<310 


0.8±0.2 


LETGS 


0±50 


<330 


0.8±0.3 


RGS 


-110 ±80 


<490 


0.9±0.1 




21.5 22 22.5 23 



Wavelength (A) 

Fig. 3. Detail of RGS 1 showing the O VII triplet as well as 
absorption lines for O VI and O V. The xabs model was used 
for fitting the warm absorber. The Galactic Ovil absorption 
line is indicated. 

are formed at the same distance from the ionizing source and 
under the same conditions, a blueshift of 620 km s _1 is incon- 
sistent with the lack of blueshift for the forbidden line. Because 
the forbidden line is strong, models predict that even tenuous 
gasses should produce a weak but detectable inter-combination 
line. One possible explanation is that the inter-combination line 
blends with an O VI absorption line at 21.79 A. 

Another emission feature for oxygen is the radiative recom- 
bination continuum (RRC) of O VII between 16.9 and 17 A (see 
Fig- El- As there is severe blending with iron absorption lines, 
we approximated the RRC simply with a delta function. We 
find a 2 a detection for this RRC and the narrowness of the 
feature (unresolved by the instrument) indicates a low temper- 
ature. The flux is 0.06 ± 0.03 ph mT 2 s _1 which is consistent 
with the 0.06 ± 0.10 ph m~ 2 s _1 value as determined from the 
LETGS spectrum (Kaastra et al. 2002a I. 

The presence of inner-shell absorption lines of oxygen to- 
gether with the higher ionized O VII and O VIII absorption lines 
allow for an accurate ionization determination, independent of 
elemental abundances. The ionization parameter determined 
for oxygen spans three orders of magnitude from log £ = —1.8 
for O in to log £ = +1.8 for O vm. 

The deepest absorption line in the spectrum is the O VIII 
Lya line. For O VIII we also detect the Ly/3, Ly7 and hy8 lines 
and some hints for the Lye line in absorption (see Fig. [2). The 
multitude and depth of these absorption lines allow thus for 
an accurate column density and outflow velocity determination 
(see Table These higher order lines are important for the 
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detection of variability in the warm absorber, as they have an 
optical depth close to unity. However, we did not detect signif- 
icant variability in the measured column densities as compared 
to the previous LETGS observation, which has rather large er- 
ror bars. Also within this observation no variability in the warm 
absorber was detected. 

Fig- shows that our model overestimates the Ovil reso- 
nance line. A probable explanation is that the emission com- 
ponent of this line partly refills the absorption line due to the 
imperfect spectral resolution. In the LETGS data (Kaastra et al. 
2002a I this line has a clear P Cygni profile, and also in the RGS 
there is excess emission at wavelengths just longer than the res- 
onance line. Both models for the warm absorber don't include 
possible P Cygni profiles. The excess emission just long-ward 
of 22.25 A is noise, as the width of the feature is narrower 
than the point spread function of the instrument. The Galactic 
absorption by Ovil, detected with 1.5 a significance is also 
indicated in Fig.|3] 

Shortward of the forbidden of O VII there is a deep O VI ab- 
sorption line, and 0.25 A longward of the forbidden line there 
is a deep O V absorption line. Both are inner-shell absorption 
lines, indicating the importance of these transitions in AGN. 
For O V two absorption lines are detected. However, the line 
at 22.33 A is redshifted by about 700 km s _1 relative to the 
weaker line at 19.92 A, and the other oxygen absorption lines. 
To check the possible uncertainty on the wavelength for the 
strongest O V line, the O V wavelength and oscillator strength 
were also calculated using the Cowan code dl981> (Raassen, 
private communication). This gives a rest wavelength of 22.38 
A, closer to the observed rest wavelength of 22.39 A, and indi- 
cates that uncertainties in the rest wavelength can be of order 
0.05 A; i.e. equivalent to the instrumental FWHM. 

4.5. Absorption by iron 

Iron is the best indicator for the ionization structure of the warm 
absorber, as it has absorption lines in the RGS band from neu- 
tral Fe I to the highly ionized Fe XXIV. For NGC 5548 we de- 
tect absorption from Fe VI to Fe XXIV (see Table [3}- The low- 
est ionization component is represented by the UTA (Sako et 
al. 2001 1 formed from inner-shell transitions of Fe VI through 
Fexi. These ions represent the lower ionized iron ions. Fig.|4] 
gives a detailed view of the spectrum between 16 and 17.5 A, 
where the UTA is located. Note that although there are some 
individual absorption features visible, the dominant effect of 
the UTA is the depression of the continuum due to an unre- 
solved blend of transitions with small cross-sections. The UTA 
clearly shows the importance of these inner-shell transitions in 
determining the continuum correctly. 

The moderately ionized component is represented by 
Fe XIII through Fe XIX, where we only find an upper limit for 
FeXH and FeXlV. The highest ionized component is repre- 
sented by absorption by Fe XIX through Fe XXIV. For Fe XIX 
through Fe XXII we detect a multitude of absorption lines be- 
tween 8 (lowest wavelength that we included in the fit) and 15 
A. For Fe XXIII and Fe XXIV only two absorption lines are de- 
tected, resulting in uncertain velocity determinations. 




x x > 

qj 0) c 



Wavelength (A) 

Fig. 4. Detail of the first order data showing the UTA between 
16 and 17.5 A. The thick line is the best fit through the data, the 
thin line indicates the same model if we set the iron abundance 
for the lowest ionization component to 0. A dead pixel occurs 
at 16.4 A. The gaps in the model at 17.05 and 17.15 A are due 
to missing data. Only the deepest features are labeled. 



It is thus clear that all three ionization components as fitted 
with xabs are amply confirmed by iron absorption. The ques- 
tion arises whether a continuous ionization structure would be 
a better fit to the data. From Fe VI through Fe XXIV we detect 
all iron ions with the exception of Fe VII, Fe XII and Fe XIV. 
Between Fe XII and Fe XV all column densities are low. Forcing 
Fe XII to Fe XVI to have a similar column densities as Fe XI or 
Fe XVII worsens the fit, increasing x„ 2 by 0.24. This indicates 
that for log £ between 1 and 1 .4 (£ between 1 and 25 ) (in 1 ~ 9 
W m) the column densities are lower than for the other ioniza- 
tion states. However, no clear depression in column densities 
is seen for the transition between the middle and highly ion- 
ized component. A more continuous ionization structure than 
presented by the xabs model for the warm absorber is not ex- 
cluded by the data. 

4.6. Absorption by carbon and nitrogen 

For carbon we detect a strong C VI Lya line (see Fig. |5}, the 
Ly/3 line (see Fig.0, Ly7 and hints of Ly<5 lines (see Fig.|6j. 
These higher order lines are important in the correct determi- 
nation of the column density as well as any possible time vari- 
ability, as their optical depth is around unity. C VI and N VII are 
produced mainly in the middle ionization component, while 
C V and N VI are produced in the middle and lowest ioniza- 
tion component. Finally, C IV and N V are produced only in 
the lowest ionization component (see Table [3] for an indica- 
tion of the ionization state). The C VI Lya line is rather deep 
and at a relatively long rest wavelength of 33.736 A, resulting 
in a very good velocity shift determination. The velocity com- 
ponent found from the X-rays is consistent with component 3 
from the UV, which has an outflow velocity of —540 km s _1 
instead of — 5 1 km s ~ 1 . There is a hint in our data that the line 
has substructure, consistent with the earlier LETGS results. 
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Fig. 5. Detail of first order RGS 1 and 2 data, indicating the 
C VI Lya absorption line as well as some C V, S XII, Ar IX, 
Arx and Arxi absorption lines. A dead pixel is seen at 34 and 
35.8 A. 




23 24 25 26 27 

Wavelength (A) 



Fig. 6. Detail of first order RGS 1 and 2 data, indicating the 
N VII Lya absorption line as well as the C VI Ly7 and Ly<5 lines. 
Absorption of S XIII, Ar XIII, and Ar XIV are labeled. Dead pix- 
els are seen at 23.4 A, 24.3 A, 25.7 A, 26.2 A, 26.6 A, and 27.6 

A. 

4. 7. Ftelativistically broadened emission lines 

After finding a model for the warm absorber that is as com- 
plete as possible, we added relativistically broadened emission 
lines of O VIII, N VII, and C VI Lya to the model. Such lines 
were observed in MCG— 6-30-15 and Mrk 766 (Branduardi- 
Raymont et al. l200Tl Sako et al. l2002l 

In NGC 5548 Kaastra et al. J2002at found evidence for 
weak broadened emission lines of O VIII Lya and N VII Lya. 
However, no indication for a relativistically broadened Cvi 
Lya line was detected. We applied the best fit parameters found 
from the LETGS model to the RGS data but found no sig- 
nificant improvement for any of the relativistically broadened 
emission lines. Leaving the parameters for the Laor line profile 
(Laor ll99"fl free, we find that a broad Gaussian would better 




28 29 30 31 32 

Wavelength (A) 



Fig. 7. Detail of first order RGS 1 and 2 data, indicating an N VI 
absorption line as well as the C VI Ly/3 line. Further absorption 
by S XII, S XIV and Ar XIII are labeled. Dead pixels are seen at 
28.1 A and 30.3 A. A CCD gap related feature is seen at 28.9 

A. 



Table 7. Values for the relativistically broadened emission lines 
in NGC 5548, MCG-6-30-15 and Mrk 766 (Sa ko et al . l200fl . 
LETGS measures are taken from Kaastra et al. (|2002a). 





RGS 


LETGS 


MCG-6-30-15 


Mrk 766 




46±35 


46(-3,+8) 


38.5±0.4 


34.3±0.74 




3.9±5 


3.9±0.6 


4.49±0.15 


3.66±0.22 


l'in C 


100±30 


<2.6 


3.21±1.2 


<2.3 


.. c 
lout 


400±90 


>10 


100(-48,+95) 


81(-20,+73) 


Cvi d 


<0.8 


<0.7 


2.0±0.2 


2.3±0.2 


Nvil d 


<0.03 


1.1 ±0.4 


4.4±0.2 


5.8±0.3 


Ovra d 


0.2±0.1 


0.6±0.2 


13.1±0.6 


11.4±0.7 



a The inclination derived for the AGN from the relativistic lines, in 
degrees. 

b The emissivity slope. 

c InGM/c 2 . 

d EW measured in A. 



describe the line/excess. However, the improvement in fit is 
still small for all three emission lines. The results are given 
in Table [7] together with the best fit values for the LETGS 
data of NGC 5548 (Kaastra et al. I2002a> . and the RGS data 
of MCG-6-30-15 and Mrk 766 (Branduardi-Raymont et al. 
2001 Sako et al. l2002l Note, that the values found for the EW 
between our RGS results and the LETGS results are consistent 
within 3 a, except for the inner radius. These relativistic lines 
should be time variable, considering that they are formed in the 
accretion disk close to the black hole. Possibly, time variability 
explains the differences between the previous LETGS results 
and ours. Given their weakness, in the further analysis and in 
all the plots we did not include these broadened emission lines. 
Fig.|H]gives the wavelength band where the relativistic O VIII 
Lya line should be situated, and a small excess is seen between 
18 and 18.4 A. 
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Wavelength (A) 

Fig. 8. Best xabs fit to the first order data, detailing the wave- 
length band around the position of the relativistic O VIII Lya 
emission line. Note that in this plot no relativistic line is fitted. 
Dead pixels are seen at 17.4 A, 18.3 A, 18.7 A, 18.8 A, 19 A, 
and 19.95 A. A CCD gap related feature can be seen between 
20.7 and 20.8 A. 



5. Discussion 

5. 1. Outflow velocity versus ionization 

No evidence is found for a correlation between the outflow ve- 
locity and the ionization parameter (see Fig. [9}. Rather the out- 
flow velocity, for first order, is consistent with —380 kms -1 
and —350 kms -1 for the slab and xabs model respectively 
over the three orders of magnitude sampled in ionization scale. 
This is consistent with the UV data (Crenshaw & Kraemer 
I1999> . and the results for the C VI Lya line in the LETGS data 
(Kaastra et al. 2002a I. Kaastra et al. ( 2002a I found a trend of 
decreasing outflow velocity for higher ionized iron and oxy- 
gen ions, although all outflow velocities were consistent with 
an outflow velocity of —340 km s _1 . This result is not repro- 
duced here. The outflow velocities observed in both data sets 
are consistent with the five UV velocity components, suggest- 
ing a similar origin of both absorbers. 



5.2. Hydrogen column density versus ionization 

The derived absorbing column density per ion, as determined 
from the slab modeling, versus the ionization parameter is plot- 
ted in Fig. Ei| In the Figure we use the ionization parameter 
corresponding to the state with maximum column density. In 
reality, however, the ion is produced over a range of ionization 
states. Fig. ^| thus assumes a continuous ionization distribu- 
tion, and is dependent on the spectral energy distribution (SED) 
in our models with XSTAR. We tested a model with and with- 
out a UV bump, and found a maximum displacement in log £ of 
0.3 for very highly ionized iron, and a maximum change in the 
derived absorbing column density of 0. 1 on the log scale used. 
However, as most ions at a particular ionization state have a 
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Fig. 9. The velocity shift as given in Table[3]versus the the ion- 
ization parameter. The vertical line represents the average out- 
flow velocity. 
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Fig. 10. The derived absorbing column density per ion as de- 
rived from the column densities given in Table [3] There is an 
increase of the derived column density with higher ionization. 



similar shift, the derived abundances and the derived absorbing 
column density are rather insensitive to the input SED. 

In general, the derived absorbing column density increases 
for higher ionization parameters. The correlation coefficient, 
for all data points, excluding upper limits, is 0.71, and the re- 
lated significance is 0.999997; for iron these values are 0.82 
and 0.9982, respectively. The derived absorbing column den- 
sity from the three component xabs model is very similar to the 
absorbing column density as derived from the slab model, and 
a similar increase for higher ionization is observed. 

In general the abundances found are consistent with so- 
lar. An exception is iron, for which both models have lowly 
ionized iron overabundant compared to lowly ionized oxygen, 
while a reverse trend occurs for highly ionized iron and oxy- 
gen. As carbon and nitrogen are consistent with oxygen, the ef- 
fect is probably due to uncertainties for the lowly ionized iron 
ions. A similar overabundance for iron in the lowest ioniza- 
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tion state was noted by Blustin et al. (2002i for NGC 3783. In 
the LETGS observation of NGC 5548 (Kaastra et al. I2002a> 
found a similar abundance pattern for iron, throughout the ion- 
ization range. The explanations listed by Kaastra et al. ( 2002a I 
for this discrepancy have been tested using this data set. First, 
a calibration effect is unlikely as RGS and LETGS give con- 
sistent results. Secondly, we determined the column density for 
lowly ionized oxygen ions from our own X-ray spectrum; and 
thus we do not need to rely on non-simultaneous UV measure- 
ments. Spectral variability is thus not a possible explanation for 
the apparent overabundance of lowly ionized iron. Thirdly, we 
tested the sensitivity of the model to different values for the line 
broadening, ensuring that saturation due to narrow line widths 
cannot be responsible for the overabundance of lowly ionized 
iron. Reducing the width of the lowly ionized absorption lines, 
has the effect that iron saturates first, worsening the overabun- 
dance already found. The maximum width of the lowly ionized 
lines is set by the width of the inner-shell oxygen lines, which 
are unresolved. 

Possible explanations for the iron overabundance are an un- 
certainty in the wavelength of the M-shell iron absorption lines; 
or if the optical depth, r > 1, an uncertainty in the radiative 
transfer models used. Another possibility is the omission of 
certain processes like re-emission and Bowen fluorescence in 
the radiative transfer models used. Finally, the iron abundance 
could be non-solar, although the physical process causing such 
an overabundance of iron in Seyfert galaxies is unknown. Non- 
solar abundances have been detected in IRAS 13349+2438 
(Sako et al. l200fl and NGC 1068 (Kinkhabwala et al. l2002l 

5.3. Outflow geometry 

The geometry and physics of the warm absorber as observed 
in some AGN is poorly understood. One of the outstanding 
questions is whether this outflow is spherical, in which case 
the ionization of the gas is mainly dependent on the distance. 
Alternatively, the outflow could be in localized streams formed 
due to instabilities in the accretion disk. The ionization param- 
eter is then mainly dependent on the density. As distance and 
density cause a degeneracy in the determination of the ioniza- 
tion parameter, the location of the outflowing X-ray and UV 
warm absorber is uncertain. Constraining the geometry of this 
outflow will help solve these outstanding problems. 

From mass conservation (A/i oss = 4 ■Kr 2 nm p vVt), the ion- 
ization parameter (£ = L/nr 2 ), the outflow velocity (v), and 
the bolometric luminosity (L), one can determine the mass loss 
rate, M\ oss , as a function of the opening angle (SI) of the out- 
flowing wind. This is summarized in eq. ([0, where m p is the 
proton mass. Assuming that the mass loss rate of the wind can- 
not be higher than the mass accretion rate onto a black hole 
and that the system is stationary, one can equate both and ob- 
tain an upper limit for the opening angle. Eq. (0 gives the mass 
accretion rate, M acc , with a mass conversion efficiency of n. 



los; 



Macc — 



C 2 T] 



(1) 



(2) 



In the case of a Schwarzschild black hole, as considered here, 
f] = 0.057. For our average outflow velocity of 380 km s _1 the 
upper limits of the opening angle vary between 3.8 x 10" 7 sr 
for log £ = -1.8 and 0.024 sr for log £ = 3, with log £ in 10~ 9 
W m. For an average ionization state of log £ = 1.7 in 10 -9 
W m, we find an upper limit to the opening angle of 0.0019 
sr. The upper limits to the opening angle indicate that the out- 
flow is mostly in a narrow stream, where the densest part of the 
stream is the narrowest and lowest ionized. For higher outflow 
velocities or a higher mass conversion efficiency, i.e. in the case 
of a Kerr black hole, the stream is even narrower. 

From Fig. ^| we determine that the derived absorbing col- 
umn density scales with the ionization parameter, Ah ~ 
a ranges from 0.25 to 0.5 depending on whether one uses the 
Ah obtained for the iron or oxygen ions at the lower ionization 
states. We define Ah = nd, where n is the density and d is 
the thickness of the outflow in the line of sight. Equating Ah to 
£ a , where £ = L/nr 2 , one can determine a relation between the 
density, thickness and distance. If one further assumes that for 
a fixed distance (r) from the ionizing source, there is a range in 
densities, resulting in a range in ionization states, one finds that 
the density is related to the thickness n ~ d~ 5 / 4 or n ~ d~ 3 / 2 , 
respectively. From both simple models we conclude that the 
outflow occurs in streams with a very small opening angle for 
the dense and lowly ionized gas, and with a larger opening an- 
gle for the less dense and higher ionized gas. 

5.4. Comparison of UV and X-ray detection of O vi 

Ovi has absorption lines in both the UV (AA1032, 1038 A) 
and X-ray band (22.01 A, 19.34 A and 21.79 A). Our X-ray 
column density of Ovi of 10 20 5 m~ 2 is nearly an order of 
magnitude larger than the UV column density (10 1962 m -2 ) 
of Brotherton et al. (2002 1. However, the outflow velocities in 
the UV and the X-rays are comparable. The difference in de- 
rived absorbing column density indicates that the UV data are 
severely saturated, leading to an underestimate in the column 
density. Arav et al. (2002 ( compare the column density of O VI 
as observed in this data set and an earlier FUSE UV measure- 
ment, and conclude that the column densities could be equal 
if velocity dependent covering is considered. However, simul- 
taneous observations are necessary to exclude the possibility 
of column density variability (see Arav et al. 2002 for more 
details). Generally, in the X-rays, we find significantly more 
low ionization gas than previously deduced from UV data. The 
discovery of O III to O V absorption lines, in addition to O VI, 
gives us for the first time an estimate for the column densities 
of these ions, which are inaccessible in the UV. 

6. Summary 

We have presented here the highest signal to noise high resolu- 
tion X-ray spectrum of NGC 5548 obtained yet. The spectrum 
shows a rich structure in narrow and broad spectral features. We 
detected a very weak O VII RRC, consistent with a low temper- 
ature. Inner-shell oxygen lines, together with the higher order 
absorption lines from O VIII Lya and C VI Lya have an optical 
depth near unity, important in the detection of spectral vari- 
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ability of the warm absorber, although no spectral variability 
was detected in this data set or from comparison with the ear- 
lier LETGS spectrum. Uncertainties in rest wavelengths for the 
inner-shell ions complicate the study of these absorption lines, 
certainly for the lowly ionized ions. Possibly this or uncertain- 
ties and the omission of certain processes in the radiation trans- 
fer models lead to an deviant iron to oxygen abundance ratio. 

The detected warm absorber spans three orders of magni- 
tude in ionization, from ions as lowly ionized as those mea- 
sured in the UV band, to nearly fully ionized iron. Only a few 
percent of the X-ray warm absorber column density is lowly 
ionized, represented by oxygen and iron inner-shell transition 
lines. The outflow velocities determined from the X-rays are all 
consistent with those measured in higher resolution UV spec- 
tra, and no correlation between the ionization state and velocity 
was detected. The similarity in outflow velocity and the detec- 
tion of O VI give strong evidence that the UV and X-ray warm 
absorber are different manifestations of the same outflowing 
wind phenomenon. However, we detect an order of magnitude 
more O VI in the X-rays than previous UV measurements. In 
general, we detect substantially more low ionized gas, than is 
deduced to be present from UV observations. Probably the UV 
band underestimates the true column densities due to satura- 
tion and velocity dependent covering factor. From simple ar- 
guments, we found that the outflow detected is consistent with 
a narrow stream, with smallest opening angles for the densest 
and lowest ionized gas. 
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